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[i] This study investigates effects of drizzle and cloud horizontal inhomogeneity on cloud 
effective radius ( r e ) retrievals from the Moderate Resolution Imaging Spectroradiometer 
(MODIS). In order to identify the relative importance of various factors, we developed a 
MODIS cloud property retrieval simulator based on the combination of large-eddy 
simulations (LES) and radiative transfer computations. The case studies based on synthetic 
LES cloud fields indicate that at high spatial resolution (~ 1 00 m) 3-D radiative transfer 
effects, such as illumination and shadowing, can induce significant differences between 
retrievals of r e based on reflectance at 2. 1 //m (r e 2 .i) and 3.7 //m (r e 3 7 ). It is also found that 
3-D effects tend to have stronger impact on r e2 .i than r e , 3.7, leading to positive difference 
between the two (Ar e-3 7 _ 2 .i) from illumination and negative Ar e , 3 . 7 _ 2 .i from shadowing. 

The cancellation of opposing 3-D effects leads to overall reasonable agreement between 
r e2 .i and r e 3 7 at high spatial resolution as far as domain averages are concerned. At 
resolutions similar to MODIS, however, r e 2 . 1 is systematically larger than r e 2 7 when 
averaged over the LES domain, with the difference exhibiting a threshold-like dependence 
on both r e2A and an index of the sub-pixel variability in reflectance (H a ), consistent with 
MODIS observations. In the LES cases studied, drizzle does not strongly impact r e 
retrievals at either wavelength. It is also found that opposing 3-D radiative transfer effects 
partly cancel each other when cloud reflectance is aggregated from high spatial resolution 
to MODIS resolution, resulting in a weaker net impact of 3-D radiative effects on r e 
retrievals. The large difference at MODIS resolution between r e 3 7 and r e2A for highly 
inhomogeneous pixels with H n > 0.4 can be largely attributed to what we refer to as the 
“plane-parallel r e bias,” which is attributable to the impact of sub-pixel level horizontal 
variability of cloud optical thickness on r e retrievals and is greater for r e2A than r e7>7 . 

These results suggest that there are substantial uncertainties attributable to 3-D radiative 
effects and plane-parallel r e bias in the MODIS r e 2 , retrievals for pixels with strong 
sub-pixel scale variability, and the H a index can be used to identify these uncertainties. 
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1. Introduction 

[ 2 ] Low-level maritime water clouds have a strong cloud 
radiative forcing [ Hartmann et al., 1992] and are thought to 
be particularly sensitive to aerosol influences owing to their 
low altitude. A key parameter in aerosol-cloud interactions is 
the cloud droplet effective radius r e , which also has a strong 
influence on cloud radiative effects. In satellite-based retrie- 
vals, r e is most frequently derived together with cloud optical 
thickness (r) from reflectance measurements at two wave- 
lengths [Nakajima and King, 1990], one usually in the visible 
or near-infrared spectral region (for example, 0.86 /im), 
where water has negligible absorption and therefore cloud 
reflection is determined mainly by r, and the other in the 
shortwave infrared (for example, 1.6 /im, 2. 1 /im, or 3.7 /im), 
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where liquid- and ice-phase water have significant absorption 
and cloud reflectance decreases with increasing r e . This so- 
called bi-spectral method has been widely adopted [ Han 
et al., 1994; Platnick and Twomey, 1994; Nakajima and 
Nakajima, 1995; Platnick et al., 2003; Roebeling et al., 
2006] and forms the basis of global surveys of particle size 
by the MODIS instruments [Platnick et al., 2003]. Note that 
in both the operational MODIS cloud product and in this 
study, r is defined with respect to the visible spectral region 
and remains almost invariant from the visible to shortwave 
infrared regions of interest in this study. 

[3] Because of their strong radiative forcing [Hartmann 
et al., 1992], even small biases in r e or r retrievals for 
low-level maritime water clouds can lead to significant 
errors in calculations of cloud radiative forcing. It is esti- 
mated that even a moderate perturbation to the r e of these 
clouds can lead to a global radiative forcing perturbation of 
around 1~2 W m -2 [Oreopoulos and Platnick, 2008]. It is 
therefore critical to identify the sources and magnitudes of 
bias in the cloud property retrievals for these clouds. 

[ 4 ] Recently, several studies [Nakajima et al., 2010b; 
Seethala and Horvath, 2010; Zhang and Platnick, 201 1] have 
shown that r c retrievals based on the 3.7 pm MODIS band 
(r e 3 7) are systematically smaller than those based on mea- 
surements in the 2. 1 pm band (r e2 .i). As shown by Zhang and 
Platnick [2011], the difference (Ar e 3 7 _ 2 .|) is strongly 
dependent on cloud regime: over the trade wind cumulus 
cloud region, where clouds are often broken, the difference can 
be as large as 10 pm. These findings indicate substantial 
uncertainties in the current satellite-based cloud r e retrievals. 

[5] Two hypotheses have been proposed to explain these 
large differences. One possibility is the presence of drizzle 
within the cloud. The weighting function for cloud reflectance 
is a function of the magnitude of absorption [e.g., Platnick, 
2000; Zhang et al., 2010], such that r e37 is weighted more 
toward cloud top than r e2A because of the stronger water 
absorption and smaller penetration depth of the 3.7 pm band. 
Drizzle in marine boundary layer clouds causes r e to increase 
toward cloud base, a vertical structure that might potentially 
produce larger values of r e2 . 1 than r e37 [e.g., Chang and Li, 
2002; Nakajima et al., 2010a, 2010b; Suzuki et al, 2010]. 

[6] The second hypothesis centers on the role of horizon- 
tal inhomogeneity, and possibly three-dimensional radiative 
transfer effects, in determining cloud reflectance [Boeke, 
2009; Hayes et al., 2010; Seethala and Horvath, 2010; 
Zhang and Platnick, 2011], Most satellite retrievals, 
including the bi-spectral method, assume that clouds are 
both horizontally and vertically homogeneous within the 
instrument field of view (typically on the order of a few 
hundred meters to a few kilometers). Unfortunately, this 
assumption can cause significant errors when the clouds 
exhibit sizable spatial variability within the sensor field of 
view [Marshak et al., 2006; Liang et al., 2009; Di Girolamo 
et al., 2010]. Zhang and Platnick [2011] noticed that the 
r e 3 7 and r e2 1 differ most in pixels with both large r e2A and 
large sub-pixel inhomogeneity, and suggested that effects 
associated with cloud horizontal inhomogeneity may play an 
important role in causing the observed Ar e31 _ 2A (see 
section 2 for details). 

[7] It must be noted that the two hypotheses should not be 
considered mutually exclusive as they may both play a role 
in the observed Ar ei37 _ 2 3 . Nevertheless, these hypotheses 


have very different implications. If Ar e 3 7 _ 2 . 1 is a result of 
drizzle, it could provide useful information, including 
remote drizzle detection [Chang and Li, 2002, 2003; 
Kokhanovsky and Rozanov, 201 1]. If Ar e 3 . 7 _ 2 .i results from 
cloud horizontal inhomogeneity, however, it reflects the 
limitations of the retrieval algorithm assumptions and might 
be used to assess retrieval quality. 

[8] Although an understanding of the relative importance of 
drizzle and cloud inhomogeneity on MODIS r e retrievals is 
desired, it is difficult to separate the influences of these factors 
using MODIS observations alone. Fortunately, high-resolution 
cloud resolving models provide a useful tool for cloud remote- 
sensing studies. A high-resolution cloud resolving model, 
coupled with a plausible treatment of cloud microphysics, can 
provide detailed information on cloud macrophysical and 
microphysical structure. Cloud fields from such a model can be 
used as input to drive 3-D radiative transfer models to simulate 
satellite observations. Then, cloud property retrievals based on 
simulated observations can be compared with the original 
simulated cloud fields to identify the influence of various fac- 
tors on passive cloud property retrievals. One of the advantages 
of such a simulator is its flexibility to accommodate various 
mechanisms at different levels of complexity. For example, the 
impact of 3-D effects can be estimated by comparing 3-D with 
1-D radiative transfer simulations. The impact of drizzle can be 
estimated by artificially removing it from the simulated cloud 
fields. In previous studies, the difference between r^ 21 than 
r e 3 7 has been assessed and analyzed based on MODIS obser- 
vations and hypotheses have been proposed to explain the 
difference. In this study, these hypotheses will be tested using 
synthetic cloud fields from a large-eddy simulation (LES) 
model. For this purpose, we have developed a MODIS cloud 
property retrieval simulator using a LES model combined with 
1-D and 3-D radiative transfer models. Below we present sev- 
eral case studies based on this simulator and discuss the 
implications of the results. 

[ 9 ] The objectives of this study are twofold. First, we will 
introduce our simulator and demonstrate its capabilities for 
simulating MODIS cloud property retrievals, including the 
subtle difference between r e2 .i and r &37 . Second, using the 
simulator, we will investigate the effects of cloud horizontal 
inhomogeneity and drizzle on satellite r e retrievals based on the 
bi-spectral method, their relative importance under various cir- 
cumstances, and the implications for introducing systematic 
bias in the observational record. In particular, we will attempt to 
address the following questions: How and to what extent do 
drizzle and cloud inhomogeneity cijfect MODIS r e2A retrievals? 
Do they’ affect the r e37 retrieval in different ways or to different 
degrees? Finally, can they lead to systematic retrieval bias? 

[ 10 ] We will first introduce in section 2 the “plane-parallel 
r e bias” using simple idealized cases, as well as summarize 
recent observational results with respect to differences between 
MODIS r e 3 7 and r e2A for maritime water clouds. We will 
then describe the MODIS simulator in section 3 and the case 
studies in section 4. Results are discussed in section 5 in the 
context of MODIS observations. 

2. Background 

[ 11 ] To provide clarity for more detailed discussions 
below, it is important to briefly introduce the concept of the 
plane-parallel r e bias and revisit some 3-D radiative transfer 
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Figure 1. Two theoretical cases to illustrate the nonlinearity effect in r e retrievals resulting from sub- 
pixel cloud inhomogeneity. Numbers on top of the Nakajima-King look-up-table (LUT) curves corre- 
spond to values of r contour lines in the LUT, and the numbers on the right of the curves correspond 
to values of r e contour lines in the LUT. 


effects known to have strong impacts on cloud property 
retrievals. A summary of the main findings from a recent 
study of MOD1S r e retrievals for maritime water clouds 
follows. 

2.1. Plane-Parallel r e Bias 

[ 12 ] We define the plane-parallel r e bias as the impact of 
small-scale variability in r on r c retrievals that use area- 
averaged reflectance. This bias is illustrated using two ide- 
alized examples in Figure 1, which shows forward calcula- 
tions of reflectance in a single non-absorbing band (0.86 / an) 
and two bands with different amounts of liquid water 
absorption (2.1 n m and 3.7 /im). In Figures la and lb, we 
assume that half of a MODIS pixel overlying a black surface 
is covered by a cloud with T\ = 2.8 and r e = 8 /im, and the 
other half is covered by a cloud with r 2 = 30.8 and r e = 8 /im. 
Figures lc and Id assume the same optical thickness but use 
r e = 18 /im. Focusing on the r retrieval, the figure illustrates 
the well-known “plane-parallel albedo bias” [Cahalan et al., 
1994]: the retrieved r based on the mean reflectance of 
inhomogeneous pixels tends to be smaller than the mean of 
the sub-pixel r . In this example the value of r retrieved from 
the averaged reflectance (10.8) is substantially smaller than 
the average value (16.8). 

[ 13 ] This problem is more acute for retrievals of r e because 
the reflectance used to infer r e also depends on r over much 


of the range of plausible values. If the reflectance at non- 
absorbing and absorbing wavelengths depended only on r 
and r c respectively, (in other words, the look-up-table is 
orthogonal) reflectance at absorbing wavelengths would be 
uniform in our example and particle size could be retrieved 
perfectly. As the figure demonstrates, however, the look-up- 
table is not orthogonal. The nonlinearity leads to a simulta- 
neous underestimate of r (i.e., plane-parallel-albedo bias) 
and overestimate of r e (i.e., plane-parallel r e bias). The area 
over which this is true is larger in the less-absorbing band, 
which explains why the size overestimate at 2.1 /im is larger 
than at 3.7 /im (r e estimates of 11.7 and 8.8 /im, respec- 
tively, in Figures la and lb). The impact becomes more 
pronounced as r e increases: in Figures lc and Id the true r e = 
1 8 /im while r e2 1 and r e 3 7 retrieved from averaged reflec- 
tances are 24 /im and 18.1 /im, respectively, resulting in a 
Ar ej 3 . 7 _ 2.1 around —6 /tm. 

[ 14 ] It is to be noted that the underlying assumption 
behind the examples in Figure 1 is that r has stronger small- 
scale horizontal variability than r e . This assumption appears 
reasonable, as r can vary over several order of magnitude, 
while r e varies mostly from a few to a few tens of microns. 
On the other hand, if r e has substantial sub-pixel variability, 
the configuration of r and r e would be different from those 
in Figure 1. The plausibility of the example shown in 
Figure 1 will be examined in the LES cases. 
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Figure 2. Dependence of Ar e>3 7 _ 2 , on cloud effective 
radius (r e2 1) and cloud horizontal inhomogeneity index 
( H a ) derived from one month of operational MOD1S obser- 
vations of warm (cloud top temperature >273 K) liquid- 
phase clouds over ocean [Zhang and Platnick, 2011]. The 
color shading, with the scale given on the right, corresponds 
to monthly mean values of Ar e 3 _ 7 _ 2 .i for clouds with r > 5 
on the space specified by r e2 i and H a . The gray lines indi- 
cate the relative frequency of each grid box, specified by 
selected combinations of r e2A and H a (unity corresponds 
to the most frequently observed combination of r e2 i and 
H a ). Thin clouds with r < 5 are excluded from the figure 
because of their substantial retrieval uncertainties. 


2.2. 3-D Radiative Transfer Effects 

[is] Cloud horizontal inhomogeneity can also induce net 
horizontal transfer of radiation that is neglected by the one- 
dimensional models on which retrievals are based. Effects of 
this 3-D radiative transfer, such as illuminating and shadow- 
ing, are known to have potential impacts on r e and r retrievals 
[e.g., Varnai and Davies, 1999; Vcirnai and Marshak, 2002; 
Kato et al, 2006; Marshak et al., 2006]. For example, the 
illuminating effect makes clouds appear brighter than expected 
under the plane -parallel cloud assumption, resulting in over- 
estimated r and underestimated r e retrievals. The shadowing 
effect makes clouds appear darker, resulting in underestimated 
r and overestimated r e [Marshak et al., 2006], 

[16] The 3-D radiative effects may affect r e2 . 1 and r e3 7 
retrievals to different extents and lead to significant differences 
between the two. In a simple step-cloud case, Zhang and 
Platnick [2011] noted that the 3-D radiative effects tend to 
have a stronger impact on r e2 1 than r e 3 7 (see their Figure 15). 
This is likely because stronger absorption in 3.7 /an leads to 
less multiple-scattering and, as a result, less photon horizontal 
transport. The implication is that the illuminating effect tends 
to result in positive Ar e 3 7 _ 2 .i because it reduces r e2 1 more 
than r e 3 7 , while the shadowing effect tends to result in nega- 
tive Ar e 3 .7-2.1 because it increases r e 2A more than r e37 . 

[17] When a cloud pixel has strong horizontal heteroge- 
neity, both 3-D radiative effects and the plane-parallel r e bias 
discussed in the previous section may collude to create either 
positive or to either positive or negative Ar e>3 . 7 _ 2 .i. How- 
ever, because the illuminating and shadowing effects natu- 
rally come in pairs and tend to cancel each other [Marshak 
et al., 2006], the plane-parallel r e bias might be expected 


to be the dominant factor when averaged over many cloud 
cells. We will return to this point in section 5.2. 

2.3. MODIS Observations 

[is] Many studies have noted that MODIS retrievals of r c2 , 
tend to be systematically larger than r e37 [Chang and Li, 2002; 
Nakajima et al., 2010b; Seethala and Horvath, 2010; Zinner 
et al., 2010] and that the difference between them, Ar e 3 7 _ 2 ls 
is a strong function of cloud regime, varying from 0 to —2 /im 
over coastal stratocumulus to as large as —5 to —10 /im in 
regions of broken cumulus [Zhang and Platnick, 2011], Zhang 
and Platnick [201 1] demonstrated that underlying the regional 
dependence was a dependence on the sub-pixel scale vari- 
ability of non-absorbing reflectance H a [Liang et al, 2009]: 

stdev[R,(0.86 /im, 250m)] 

mean[R,(0.86 /im, 250m)] ’ 

where stdev[ft,(0.86 /an, 250 m)] and mean[/?,(0.86 /an, 
250 in)] indicate the standard deviation and mean of the 
measured reflectances, respectively, for the principle sixteen 
250 m-resolution sub-pixels within the 1 km MODIS retrieval 
footprint. Thus, H a has a spatial resolution (i.e., 1 km) con- 
sistent with the cloud property retrieval and increases with 
pixel inhomogeneity. The dependence of pixel-level 
Ar e 3 7— 2 1 on r e 2 1 and H a in one month of MODIS obser- 
vations of warm liquid phase clouds (with cloud top temper- 
ature >273 K) over ocean is illustrated in Figure 2. The most 
striking feature is the area of large differences (Ar e 3 7 _ 2 i ~ 
— 10 f im) in the upper right comer of Figure 2. The large dif- 
ferences are associated with large (>25 /im) values of r e2A and 
large sub-pixel inhomogeneity (H a ~ 1). Ar e 3 7 _ 2 | remains 
stable when II „ is smaller than about 0.4, but becomes 
increasingly negative with increasing H a when H a >0.4. For a 
given H a , Ar eA7 _ 2A generally increases with r e2A . 

[19] As mentioned earlier, some studies interpret negative 
values of Ar e 3 7 _ 2 | as evidence of drizzle in water clouds, 
such that the 2.1 /xm band (with weak water absorption) 
penetrates deeper into the cloud than the 3.7 /a n band and so 
can be more strongly affected by large drizzle drops in the 
lower part of the cloud [Chang and Li, 2002; Nakajima 
et al., 2010a]. Others argue that the difference between 
r e 2A and r e 3 7 is mainly a result of the plane-parallel r e bias 
illustrated in Figure 1. Both hypotheses would find support 
in Figure 2. The fact that the strongest Ar e 3 7 _ 2 | is found 
where r e 2A is larger than 25 /an is consistent with the drizzle 
hypothesis because drizzling clouds tend to have larger r e 
than non-drizzling clouds. On the other hand, large values of 
the sub-pixel inhomogeneity index II„ in the same region 
would imply plane -parallel r e bias in the retrieval. 

3. Simulating MODIS Retrievals 

[20] Since the two hypotheses cannot be distinguished 
using the observations alone, we turn to retrievals from 
synthetic clouds in which the radiative transfer effects can be 
carefully controlled. Cloud structure is obtained from simu- 
lations by the DHARMA LES model [Stevens et al., 2002; 
Ackerman et al., 2004], Here, DHARMA uses a single- 
moment bin microphysics scheme to resolve size distribu- 
tions of aerosols and activated water drops each into 25 size 
bins, spanning particle radii of 0.01-2.5 and 1-250 /an, 
respectively. The size distributions of activated droplets on 
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Table 1. Summary of Model Setup and Cloud Properties of the Three LES Cases Used for This Study 


Case Name 

Domain Size 

Grid Mesh 

Cloud 

Type 

Total Particle 
Cone. (#/cm 3 ) 

Cloud 

Coverage 3 

Domain 
Average r 

Domain 
Averaged 
LWP (g/m 2 ) 

Domain Averaged 
Precipitation Rate 
at LCL b (mm/d) 

ATEX Clean 

9.6 km x 9.6 km 

Ax = Ay = 100 m Az = 40 m 

Cu 

40 

72% 

5.4 

70 

0.26 

ATEX Polluted 

9.6 km x 9.6 km 

Ax = Ay = 100 m Az = 40 m 

Cu 

600 

78% 

8.6 

41 

0.00 

DYCOMS-II (RF02) 

6.4 km x 6.4 km 

Ax = Ay = 50 m Az varies 

Sc 

75“ 

100% 

14.5 

114 

0.07 


a Cloud coverage is defined as fraction of LES columns with 0.86 /jm cloud reflectance larger than 0.02. 
b Domain-average lifting condensation level (LCL) used as representative of cloud base. 

“Aerosol size distribution is himodal for this case; given is the concentration in the accumulation mode. 


the DHARMA grid is used to derive cloud scattering proper- 
ties at MODIS cloud bands, which in turn are used to drive 
radiative transfer simulations. The bin microphysics approach 
obviates the requirement of a priori assumptions about cloud 
drop size distributions. Fields from the LES model described 
by Feingold et al. [1996] and Xue et al. [2008] have also been 
examined, but are not shown here because the results are 
comparable to our analysis of the DHARMA cases. 

[ 21 ] Radiative transfer is computed using the 13RC com- 
munity Monte-Carlo model [Cahalan et al., 2005; Pincus 
and Evans, 2009] for 3-D simulations and the DISORT 
model [ Stamnes et al., 1988] for 1-D simulations. The single- 
scattering properties of cloud droplets (extinction efficiency 
Q ei , single-scattering albedo tu,- and phase function P,) in each 
size bin, i, are computed using the Mie code of Wiscombe 
[1979], following the steps described by Platnick and Valero 
[1995], for each MODIS band used. In the DISORT model, 
the scattering properties of a cloudy grid cell are computed 
off-line by averaging Q ei , tu; and I 1 , over the cloud drop size 
distribution (DSD) from the LES. In the I3RC model, the 
DSD-averaged scattering properties are constructed online 
(instead of pre -computed off-line as in the 1-D simulation), 
from Q e i , uji, and P, based on the cloud DSD. The dynamical 
construction is implemented through a Monte-Carlo sampling 
of the relative cumulative extinction function F defined as: 

i 

Qejfj n J 

F; = ^ ,i= I...N. (1) 

Qejfj n J 

j = 1 

[ 22 ] For a given MODIS band, the shape of F is deter- 
mined by the cloud DSD. In a scattering event, a random 
number (, uniformly distributed between 0 and 1, is gener- 
ated. Cloud droplets in the zth size bin with F, < ( < F l+1 are 
then chosen to interact with (i.e., scatter or absorb) the 
photon. This on-online construction scheme preserves the 
detail of the cloud microphysical properties from the LES in 
radiative transfer simulation. 

4. Case Studies 

4.1. Cloud Fields From Large-Eddy Simulations 

[ 23 ] Using the MODIS simulator described above, we 
have investigated the MODIS r e retrievals for three large- 
eddy simulations of marine water clouds. The first (referred 
to as “ATEX clean” hereafter) and second (“ATEX pol- 
luted”) cases are based on an idealized case study [Stevens 
et al., 2001] from the Atlantic Trade Wind Experiment 


(ATEX), with different initial aerosol loadings. A diagnostic 
approach [Clark, 1974] is used for the aerosol, in which the 
total number concentration of unactivated aerosol particles 
plus activated droplets is fixed at values of 40 and 600 cm -3 
for the clean and polluted cases, respectively, resulting in 
average cloud droplet concentrations (weighted by liquid 
water mixing ratio) of about 30 and 300 cm“ 3 . The total 
number and size distribution of activated droplets vary in 
each grid cell over the course of the simulation. For the 8-h 
ATEX simulations the domain size is 9.6 km x 9.6 km x 
3 km, with a uniform grid mesh of Ax = Ay = 100 m and 
Az = 40 m. (Further details of the model setup for the ATEX 
cases are provided by Fridlind and Ackerman [2011].) A 
snapshot of the cloud field is taken every half hour over the 
last 4 h. Therefore, each case contains eight 3-D cloud 
scenes from the large-eddy simulation, equivalent to about 
800 1 km x 1 km pixels. The third case (referred to as the 
DYCOMS-1I case) is based on an idealization of conditions 
observed during the second research flight (RF02) of the 
Second Dynamics and Chemistry of Marine Stratocumulus 
project (DYCOMS-II). For this case, the domain size is 
6.4 km x 6.4 km x 1.5 km, and the grid spacing is Ax = 
Ay = 50 m in the horizontal and stretched vertically, with 
a spacing minimum of 5 m near the surface and inversion. 
The average cloud droplet number concentration in the 
DYCOMS-II case is about 60 cm~ 3 . A snapshot is taken 
every hour over the last 4 h of the 6-h simulation, yielding 
about 200 1 km x 1 km pixels. (Further details of the model 
setup for the DYCOMS-II case are provided by Ackerman 
et al. [2009].) Table 1 provides a summary of model setup 
and cloud properties for the three cases. 

[ 24 ] Figure 3 provides a planar view of the cloud r derived 
from LES for these three cases. The DYCOMS-II case is 
almost overcast, while the cloud fraction, defined as fraction 
of columns with r > 0.1, in ATEX is around 0.7. Enough 
drizzle develops in the ATEX clean case that drizzle (i.e. ,r> 
30 //,m) accounts for a significant fraction of the total optical 
thickness r, but is negligible in the other cloud fields. It is 
important to note that the cloud r in Figure 3 is derived from 
the droplet number concentrations from the LES. It can be 
different from the r retrieved from cloud reflectance. The 
three sets of cloud fields also differ significantly in terms of 
r e and cloud inhomogeneity, which we will exploit to help 
understand the dependence of Ar ei3 . 7 _ 2 .i on cloud regimes. 

[ 25 ] Figure 4 provides an example cloud microphysics 
simulation from the LES model. The cloud droplet size 
distributions at three different vertical levels of a selected 
LES column from the ATEX-clean case are plotted. The 
location of this column is indicated in Figure 3 (left) by a red 
star. This particular column is drizzling as indicated by the 
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ATEX Clean 



Figure 3. Plan view of the cloud t of (left) the ATEX clean case at 6 h of simulation time, with the red 
star indicating the location of the column shown in Figure 4, (middle) the ATEX polluted case at 6 h, and 
(right) the DYCOMS-1I case at 4 h. The red contour indicates where drizzle drops with r > 30 //m contrib- 
ute more than 10% of the cloud r. 


development of the drizzle mode in the size distribution 
from cloud top toward cloud base. Figure 4 also demon- 
strates the level of complexity of the LES cloud fields used 
in this study. In many previous studies based on LES cloud 
fields, cloud microphysical properties are often greatly 
simplified. For example, in Kato et al. [2006] and Marshak 
et al. [2006] a constant value of r e is assumed for all 
cloudy cells. In contrast, in this study every detail of the bin 
microphysics from the LES model is preserved in the radi- 
ative transfer simulations as described in section 3. 

4.2. Radiative Transfer Calculations 

[ 26 ] Both 1-D and 3-D radiative transfer calculations are 
performed for the 0.86 pm, 2.1 pm and 3.7 pm MODIS 
bands. Because most of the low-level maritime warm clouds 
of interest to this study are located in tropical or sub-tropical 
regions, reflectance is calculated at two solar zenith angles 
(SZA) (20° and 50°) and a single value of relative azimuth 



Figure 4. Cloud droplet size distributions at three different 
vertical levels of a LES column (red star in Figure 3) from 
the ATEX clean case. The three levels roughly correspond 
to cloud top (altitude z = 1.5 km), middle of the cloud (z = 
1.0 km), and cloud base (z = 0.7 km), respectively. 


(30°) to simulate the high-sun and low-sun conditions fre- 
quently encountered by MODIS over the tropical or sub- 
tropical regions [Seethala and Horvath, 2010]. Reflectances 
under each SZA are computed at viewing zenith angles (10° 
intervals from —50° to 50°) that mimic the scan range of the 
MODIS instrument. Without explicit indication, retrieval 
results at different solar and viewing angles are simply 
combined together in the analysis for more robust statistics; 
the dependence of results on solar and viewing angles will 
be investigated later in section 5.2. For simplicity, the sur- 
face is assumed to be black and atmospheric absorption is 
not considered in any of the simulations. 

[ 27 ] In practice the 3.7 pm band observation contains both 
solar reflection and thennal emission components, the latter 
having to be accounted for during retrieval. We initially 
included emission in our radiative transfer calculations; 
however, sensitivity studies (not shown) indicate that 
including the thermal correction makes no significant dif- 
ference in the simulation, especially for pixels with r > 5. 
For this reason, we will consider only the solar reflection 
component of the 3.7 pm observation. 

[ 28 ] We perform retrievals using reflectance at high 
(native LES) and MODIS-like resolutions; for the latter we 
use area-averaged reflectance as input to the retrieval. (We 
have adopted 800 m resolution instead of the nominal 1 km 
because our model domains are more neatly divisible by 
800 m. Sensitivity studies [not shown] indicate the differ- 
ence between 800 m and 1 km retrievals is negligible.) 
Given two methods for computing radiative transfer (1-D 
and 3-D), this produces four sets of retrievals, each of which 
is affected by different factors (Table 2). High-resolution 1-D 


Table 2. Summary of How Different Effects Influence the Retrievals 
Based on the Different Combinations of Spatial Resolution (LES 
Native or MODIS) and Radiative Transfer Schemes (1-D or 3-D) 



3-D 

Radiative 

Effects 

Cloud Vertical 
Structure and 
Drizzle Effects 

(Plane-Parallel 
r e bias) 

1-D LES resolution 

No 

Yes 

No 

1-D MODIS resolution 

No 

Yes 

Yes 

3-D LES resolution 

Yes 

Yes 

No 

3-D MODIS resolution 

Yes 

Yes 

Yes 
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Figure 5. Joint probability distribution of r e ^j and r c2 i 
at the LES resolution for the ATEX clean case: (a) r e 3 7 
versus r e2A based on 1-D radiative transfer simulation, and 
(b) Ar e> 3 . 7 _ 2 .i versus r e 2 . i from 1-D radiative transfer simu- 
lations. (c, d) Same as Figures 5a and 5b but are based on 3-D 
radiative transfer simulations. Dotted lines indicate 1 : 1 rela- 
tion between r e ^j and r^ 2A . 


simulations, for example, are free of 3-D radiative effects and 
plane-parallel r e bias, but are affected by the cloud vertical 
structure, including drizzle, so that differences between r e2 .\ 
and /' e 3 7 retrievals can be attributed to cloud physics. Retrie- 
vals simulated at 800 m resolution based on 3-D radiative 
transfer are affected by all potential factors and resemble real 
MOD1S observations. It is worth pointing out that, no matter 
which method (i.e., 1-D or 3-D) or resolution (LES native or 
MODIS) are used in the forward radiative transfer simulations, 
to be consistent with the operational MODIS product the 
retrievals are based on the look-up-tables under plane-parallel 
cloud assumption. A simple threshold method based the cloud 
reflectance in the 0.86 pm band (i.e., R(0.86 pm) > 2%) is used 
to mask the cloud pixels after radiative transfer simulation. 
Only pixels with fractional cloud cover exceeding 0.95 are 
used in the analysis. 

[ 29 ] Before proceeding to retrieval results, it is important 
to point out that although a “reference r” is desirable for the 
comparison between r e 3 7 and r e2A , there are many ways 
(and therefore no unambiguous way) to define a reference r e 
from the LES field. To appreciate this, one can consider the 
LES column in Figure 4. A reference r e can be defined as 
r e = 3LWP / (2 rp w ), where LWP is the cloud liquid water 
path of the column and p w is the density of water. Such a 
reference r e ensures the correct LWP when r is known. 
However, this definition does not account for the fact that 
the MODIS 2.1 pm and 3.7 / / in bands are only sensitive to 
the upper part of the cloud because of their limited 


penetration depth [Platnick, 2000]. One way to account for 
this sensitivity is to first derive a vertical weighting function 
that incoiporates both the underlying cloud microphysical 
structure and the sensitivity of each MODIS band to that 
structure. Then a reference r e consistent with the MODIS r e 
retrieval mechanism can be derived, based on this weighting 
function, as described by Platnick [2000]. But such a refer- 
ence r e is equivalent to the band-specific r e retrieved from 1-D 
radiative transfer simulations at LES resolution. Therefore, we 
simply use the r e retrieval results from 1-D radiative transfer 
simulations at LES resolution as the reference r e . 

4.3. Results From the ATEX Clean Case 

[ 30 ] Figure 5 shows comparisons of 2 .i and r 6 j 3.7 retrie- 
vals for the ATEX clean case at the LES resolution based on 
1-D (Figures 5a and 5b) and 3-D (Figures 5c and 5d) radi- 
ative transfer simulations. To be consistent with the obser- 
vational results in Figure 2, only those pixels with cloud 
optical depth r > 5 are included in Figure 5. 

[ 31 ] There is relatively little variability in reference r e in 
this simulation: most values lie between 13 pm and 22 pm, 
with very little difference between estimates from 2.1 pm 
reflectance and those from 3.7 pm reflectance (Figures 5a 
and 5b). The close agreement between r c37 and r c2 1 in 
Figure 5a suggests that drizzle or cloud vertical structure has 
either little impact, or affects reflectance at both wavelengths 
roughly equally. We return to this point in section 5.1. 

[ 32 ] In 3-D simulations retrieved drop sizes agree well in 
the range between 1 5 pm to 20 pm, but r e 3 7 overestimates 
r e2A for r e 2 1 < 15 pm and underestimates r e2 t for r e,2.1 > 
20 pm. Most reference values of both r e 3 7 and r e2A are 
between 1 5 pm and 20 pm, so these out-of-range retrievals 
(i-e., r e , 2 .i < 15 and r e,2.l > 20 pm) are likely the result of 3-D 
effects. The results are also consistent with the idea that 
retrievals in the more-strongly absorbing 3.7 pm band are 
less affected by 3-D effects than the 2. 1 pm band. Under this 
hypothesis, for example, when 3-D effects lead to an over- 
estimate of true r e , the overestimate in r e2A should be larger 
than in r e3 7j yielding negative Ar e 3 7 _ 2 A . This hypothesis 
will be further explored in section 5.2. 

[ 33 ] Figure 5 indicates that, though 3-D radiative transfer 
increases the scatter, retrievals made from reflectance com- 
puted at high resolution are broadly consistent across 
wavelengths. Results at the MODIS resolution, shown in 
Figure 6, tell a different story: r e2il retrievals in both 1-D 
and 3-D simulations are biased significantly smaller com- 
pared with the r e 2 1 retrievals. Note also that some very large 
(>20 pm) r c 2 , retrievals emerge in Figure 6a; Ar c37 _ 2 , 
can be up to —10 pm in these pixels. The 3-D retrieval 
results in Figures 6c and 6d are quite similar to the 1-D 
results, suggesting the same mechanism is operative in both 
1-D and 3-D retrievals. 

[ 34 ] Our results are also consistent with the MODIS 
observations shown in Figure 2. Figure 7 shows composite 
plots of MODIS resolution Ar e37 _ 2 A as a function of r e 2A 
and PI u for both 1-D and 3-D simulations of the ATEX clean 
case. In the 1-D retrieval (Figure 7a), r e 2A is mostly between 
15 pm and 20 pm and A r e , 3 . 7 - 2 .i remains close to zero until 
H a reaches 0.5, after which r c 2 , rapidly jumps to much 
larger values and Ar e 3 7 _ 2 1 drops dramatically from near 
zero to negative values as large as —10 pm. Consequently, 
the largest negative values of Ar e 3 7 _ 2A are found at the 
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Figure 6. Same as Figure 5 but at MODIS resolution. 


upper right comer of Figure 7a, where both r e2] and H a 
are large. The 3-D retrievals (Figure 7b) behave much like 
the 1-D results (Figure 7a), though with more spread and 
with a less defined threshold. Unique to the 3-D simulations 
is the transition of Ar e3 . 7 _ 2 .i from red (positive) to green 


(negative) as r &2 l increases from 12 /im to more than 20 fj,m. 
Positive values of Ar e3 j- 2A correspond to points above the 
zero line in Figure 6d and above zero in Figure 5d. As dis- 
cussed above, these points are likely attributable to 3-D radi- 
ative effects, such as photon horizontal transport, illuminating 
and shadowing. 

4.4. Results From the ATEX Polluted 
and DYCOMS-11 Cases 

[ 35 ] Results from the ATEX polluted and the DY COMS-II 
cases are consistent with these findings (Figure 8). By design, 
drop sizes are smaller in the ATEX polluted case (Figures 8a 
and 8b) than in the ATEX clean examples, with reference r e 
values generally below 10 /rm. Most pixels in the ATEX 
polluted case have small-to-moderate II, T < 0.5 and stable 
values of r e2 1 and Ar e>3 . 7 _ 2 .i, with larger r e2 1 and large 
negative A r e 3 7 _ 2 1 for the small population of pixels with 
Ha > 0.5. 

[36] The DYCOMS-II case (Figures 8c and 8d) is almost 
overcast and the clouds are more homogenous, so that H a is 
small everywhere in the domain and neither the threshold 
behavior of r e21 nor very large negative A r^ 37 _ 2l values 
(i.e., Ar e 3 7 _ 2 j < —5 /tm) are evident. The effects of 3-D 
radiative transfer are similar to those demonstrated for the 
ATEX clean case. In particular, 3-D effects have a greater 
impact on r e 2 l than on r e3 7 so that anomalously small values 
of r & 2 .1 (here, less than 14 ;im) are associated with positive 
Ar e3 . 7 _ 2 .i, while pixels with abnormally large r e ,2.1 > 18 /im 
tend to have negative Ar e>3 . 7 _ 2 A . 

[ 37 ] Figures 7 and 8 demonstrate that differences between 
r e 3 7 and r e 2A retrievals at MODIS resolution in our simu- 
lations are largely controlled by the sub-pixel level cloud 
inhomogeneity: A r e , 3 . 7 - 2.1 becomes increasingly negative as 
H a increases, and the differences are associated with large 



Figure 7. Ar e 3 7 _ 2 ! at MODIS resolution for the ATEX clean case based on (a) 1-D and (b) 3-D radi- 
ative transfer simulations plotted against the sub-pixel inhomogeneity ( H a ) and r e 2 . 1 retrieval. Each dot in 
the plot indicates a MODIS resolution pixel with color indicating the value of Ar e37 _ 2A . 
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Figure 8. Same as Figure 7 but for (a, b) the ATEX polluted case and (c, d) the DYCOMS-II case. 


H a ~ 1, and large r c2 . i retrievals. These results agree well 
with the MODIS observations shown in Figure 2. 

5. Disentangling the Impacts of Drizzle 
and Cloud Inhomogeneity 

[ 38 ] In this section we investigate the effect of drizzle and 
sub-pixel level cloud inhomogeneity on MODIS r e retrie- 
vals, focusing on pixels with the largest magnitude of 
Ar e 3 7 _ 2 j to investigate causal links among r e21 , H a and 
Ai' e , 3 . 7 _ 2 . i - We use the ATEX clean case because of its 


significant drizzle and general consistency with the other 
two cases. 

5.1. Influence of Drizzle 

[ 39 ] We define drizzle drops as those droplets with r > 30 //m, 
as this threshold generally separates the droplet size distri- 
bution into drizzle and cloud modes (see Figure 4). We 
investigate the impact of these drops on our retrievals by 
perfonning calculations in which these drops are removed. 
Comparisons of LES resolution t, r c2] and r e3 7 retrievals 
based on 3-D reflectances between results that include or 
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Figure 9. Comparison of (a) r, (b) r e 2 .i, and (c) r e 3 7 between simulations at MOD1S resolution with 
and without drizzle (defined as drops with radius greater than 30 /urn) based on 3-D radiative transfer, 
(d, e) Same as Figure 6d and Figure 7b, respectively, except that plots in this figure are based on simula- 
tions without drizzle. 


exclude drizzle drops are shown in Figures 9a, 9b and 9c, 
respectively. The MODIS resolution retrievals after the 
removal of drizzle are shown in Figures 9d and 9e in the same 
manner as in Figures 6d and 7b. The removal of large drizzle 
drops in the radiative transfer simulation results in slightly 
smaller cloud optical thickness (Figure 9a) and effective 
radius retrievals (Figures 9b and 9c), as might be expected, 
but the differences between simulations with and without 
drizzle are very small. Indeed, at the MODIS resolution 
(Figures 9d and 9e), the retrievals after the removal of drizzle 
look almost identical to those with drizzle (Figures 6d and 7b). 
Thus drizzle drops with r > 30 /jm have a very minor impact 
on the r e2 .i and r retrievals in this particular example, so 
drizzle cannot be the primary reason for bias between r c 2 | and 
r e 3 7. Similar findings are reported by Zinner et al. [2010] and 
Painemal and Zuidema [2011]. 

[40] We emphasize, however, that the above results are 
based on a particular LES case, which is not representative 
of all MODIS observations. For example, a significant 
fraction of MODIS pixels in Figure 2 have relatively small 
H a ~ 0.1 and large Ar e 3 7 _ 2 1 ~ —5 //m. However, none of 
our LES cases yields any retrieval in this region. The impact 
of drizzle for the pixels in this region therefore remains 
unclear and needs to be elucidated in future research. 

5.2. Influences of Cloud Horizontal Inhomogeneity 

[41] We now turn our attention to the plane-parallel r e bias 
and 3-D radiative transfer effects. Using the ATEX-clean 


case we will show that both of these effects exist, and will 
demonstrate their influence the r e retrieval. We will then 
analyze their relative role in causing the Ar e 3 7 _ 2 1. 

[42] The idealized cases in section 2.1 illustrated how 
cloud horizontal inhomogeneity can result in a plane-parallel 
r e bias and significant difference between r e2 1 and r e3 7 . This 
bias is examined here in a more realistic setting using a pixel 
from the ATEX clean case. Figure 10 shows the cloud 
extinction coefficient and cloud reflectance along the cross 
section at SZA of 20° and viewing zenith angle of 0°. The 
vertical lines indicate the location of the MODIS resolution 
pixel analyzed in Figure 1 1 . The 1-D and 3-D simulations of 
cloud reflectance are in general agreement (i.e., net hori- 
zontal transport is small) in all but a few regions in 
Figure 10b in which K(().H6 //m) based on 3-D simulation is 
substantially smaller than that based on 1-D simulation. The 
selected MODIS resolution pixel (i.e., 800 m) consists of 64 
LES resolution sub-pixels (i.e., 100 m), the reflectances of 
which are mapped (blue asterisks) onto the retrieval LUTs in 
Figure 11. The values are clustered around the r e = 19 fi m 
contour line of the LUT. The red diamonds in the figure 
indicate the locations of the MODIS resolution cloud 
reflectances, which are simply the linear averages of the blue 
asterisks. As in the idealized case shown in Figure 1, the 
MODIS resolution cloud reflectances fall below the envelope 
of the sub-pixel values in both 1-D (Figures 1 la and 1 lb) and 
3-D (Figures 11c and 1 Id) simulations. As a result, retrievals 
of r e ,i. 1 at MODIS resolution are around 25 //m, substantially 
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Figure 10. (a) The cross section of cloud extinction coeffi- 
cient (/3) along y = 2 km in Figure 3 (left). Cloud bi-directional 
reflectance along the cross section is shown for the (b) 0.86 ^im, 
(c) 2.1 /im, and (d) 3.7 ;im MOD1S bands simulated using 
1-D (blue) and 3-D (red) radiative transfer models. The ver- 
tical black lines indicate the location of the selected pixel. 


larger than the sub-pixel mean value of 19 ^im, regardless of 
whether 1-D (Figure 11a) or 3-D (Figure 11c) radiative 
transfer is used. Retrievals using 3.7 /.tm reflectance are less 
strongly affected: MODIS resolution r e 3 7 retrievals (20.5 //m 
and 21 //m, for 1-D and 3-D radiative transfer simulations, 
respectively) are much closer to the sub-pixel mean. 

[ 43 ] The example in Figure 11 shows the existence of 
plane -parallel r e bias in our simulation, but is it typical? To 
investigate this question, we first constructed a sub-pixel 
variability index for both r and r e based on the retrievals 
from 3-D radiative transfer simulations at LES resolution. 
The index is defined as the ratio of standard deviation to 
mean of sub-pixel values. Within the context of Figure 1 1, a 
large sub-pixel variability index indicates that the sub-pixel 
level points are more scattered in the Nakajima-King LUT, 
whereas small sub-pixel variability index indicates that the 
sub-pixel level points are clustered around certain constant 
contour lines of r or r e . Figure 12b shows the histograms of 
the sub-pixel variabilities of r and r e for the ATEX-clean 
case (black line). Evidently, r has much stronger sub-pixel 
variability than r e . This suggests that the example in 


Figure 1 1 is a typical case and the r e heterogeneity factor is 
an important factor in causing the Ar e37 _ 2A . It is worth 
mentioning that we have also investigated the sub-pixel 
variability of r and r e in real MODIS observation. We 
developed a research level algorithm to retrieve r or r e at 
500 m using 500 m radiance observations from the MODIS 
0.86 /jm and 2.1 fim bands. This algorithm is applied to a 
MODIS-Terra granule shown in Figure 12a. The sub-pixel 
variabilities of r or r e derived from the 500 m MODIS 
retrievals (red lines) for this granule are shown in Figure 12b. 
The good agreement between MODIS observation and LES 
results confirms that r tends to vary more strongly than r e 
within a 1 km MODIS pixel. 

[ 44 ] We now turn to the 3-D radiative transfer effects. First, 
the impacts of 3-D radiative transfer effects on r e 2 .\ and r e 3 7 
are investigated in Figure 13. In the figure we use the dif- 
ference between the r retrieval based on 3-D radiative 
transfer simulation (referred to as “3-D r ”) and that based on 
1-D radiative transfer simulation (referred to as “1-D r ”) as 
an index of the 3-D radiative transfer effects. By definition, 
3-D r > 1-D t for the illuminating effect, and 3-D r < 1-D r 
for shadowing. As expected, the r e retrievals based on 3-D 
radiative transfer simulation (“3-D r e ”) appear smaller than 
those based on 1-D radiative transfer simulation (“1-D r e ”) 
for the illuminating effect (i.e., 3-D r e — 1-D r e < 0), and vice 
versa for shadowing (i.e., 3-D r e — 1-D r e > 0). What is 
interesting in Figure 13 is that the impact of 3-D radiative 
effects on r e 21 is generally greater than for r e 3 7 . This is 
likely attributable to the stronger absorption in the 3.7 ^m 
band that acts to reduce the horizontal transport of photons. 
The impact of 3-D radiative transfer effects on Ar ej3 . 7 _ 2 .i is 
investigated in Figure 14 as a function of SZA and retrieval 
resolution. At SZA = 50°, illuminating and shadowing 
effects are clearly seen in both LES (Figure 14a) and MODIS 
(Figure 14b) resolution retrievals. Interestingly, based on the 
3-D radiative transfer simulations, Ar e 3 7 _ 2 1 > 0 over the 
illuminated side, and Ar e 3 7 _ 2 .i < 0 over the shadowing side. 
These are consistent with the result in Figure 13 that illumi- 
nating and shadowing effects have stronger impacts on r e2 \ 
than on r e 3 7 . For example, over the illuminating side the 
apparent brightening decreases the r e retrieval, and this 
impact is stronger for r e2 l than r e37 , leading to positive 
Ar e 3 7 _ 2 |. This also explains the origin of the positive 
Ar ej3 7 _ 2 1 points in Figures 7 and 8. At SZA = 20°, a sig- 
nificant portion of 3-D r retrievals are smaller than their 1-D 
counterparts. This is likely attributable to horizontal photon 
transport. Unlike illuminating and shadowing effects, which 
are a result of geometrical cloud top height variation, the 
horizontal photon transport effect is a result of horizontal r 
variation. The absorption in the 2.1 /i m and 3.7 /im bands 
reduces horizontal photon transport, which explains why the 
absolute values of Ar ej3 . 7 _ 2 .i in Figure 14c are generally 
smaller than those of the low sun case in Figure 14a. 

[ 45 ] Another interesting point to note in Figure 14 is that, 
at LES resolution (i.e., Figures 14a and 14c), A r e ^ 37 _ 2] 
shows no obvious positive or negative bias. However, at 
MODIS resolution, Ar e 3 7 _ 2 ] seems biased more toward 
negative values, which is consistent with the results in 
Figures 5 and 6. This shift of A r e37 _ 2A is more clearly seen 
in Figure 15, which shows Ar ei3 . 7 _ 2 .i averaged over all 
cloudy pixels with r > 5 in the domain for different solar and 
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Figure 11. Cloud reflectances of the selected pixel in Figure 10 plotted in the Nakajima-King LUT. The 
blue asterisks indicate the reflectance simulated at the LES resolution and the red diamond indicates the 
MODIS resolution radiance calculated as the arithmetic average of the blue dots. The solar and viewing 
zenith angles for this plot are 20° and 0°, respectively, in this figure. 


zenith angles, as well as different retrieval resolutions. In 
both low sun (SZA = 50°) and high sun (SZA = 20°) 
simulations, and for all viewing zenith angles, the domain- 
averaged MODIS resolution Ar e 3 7 _ 2 . i (red lines) based on 
3-D radiative transfer simulations is systematically smaller 
than its LES resolution counterpart (black lines). This indi- 
cates that the shift of Ar e 3 7 _ 2 .i from close-to-zero values to 
more negative values when reflectance is averaged from 
high resolution to MODIS resolution is a robust result, only 
weakly affected by solar and viewing zenith angles. 


[46] The results in Figures 14 and 15 seem to suggest that 
although the 3-D radiative transfer effect has a strong impact 
on r e retrievals, it tends to result in random errors rather than 
systematic bias. Therefore it is reasonable to hypothesize 
that the systematic shift of Ar e 3 7 _ 2 | seen in Figures 6 and 
15 is mainly attributable to the plane-parallel r e bias. This 
hypothesis is further investigated in Figures 16 and 17. 
Figure 16 shows a comparison of the MODIS resolution 
retrievals with the mean of the LES resolution retrievals 
(referred to as the sub-pixel mean hereafter). In the case of 
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Figure 12. A comparison of sub-pixel variability of r and r e between results from MOD1S 500 m retrie- 
vals and those based on LES. (a) The RGB image of a Terra-MODIS granule (collected at 20:55 UTC on 
April 2nd, 2005) selected for comparison. A research level algorithm is developed to retrieve r and r e from 
500 m resolution radiance from MODIS band 2 (0.86 /tm) and band 6 (2.1 //m) for this selected granule. 
Sub-pixel variability of r and r e are derived for 1 km resolution operational MODIS pixels (i.e., each 1 km 
pixel contains four 500 m sub-pixels), (b) The MODIS 500 m retrieval results (red) in comparison with 
those from the ATEX clean case (black). 


the r retrieval, the plane-parallel bias is evident: values 
retrieved from the averaged radiances (i.e., MODIS resolu- 
tion retrieval) are significantly smaller than the sub-pixel 
mean values. However, in the case of r e 2 .\ retrievals, the 
reverse is true: MODIS resolution retrievals tend to be larger 
than the sub-pixel mean values. This difference indicates 
that in comparison with other factors, such drizzle and 3-D 
radiative transfer, the plane-parallel r e bias is the dominant 
source of errors in the r e 2 .\ retrievals at MODIS resolution. 
The relationship between the plane-parallel r e bias and 
Ar e> 3 . 7 — 2 .i is shown more clearly in Figure 17. In both 1-D 
and 3-D simulation, the magnitude of Ar e3 7 _ 2 \ is, in gen- 
eral, positively correlated with the strength of the plane- 
parallel r e bias because the plane-parallel r e bias has a 
stronger impact on the r e2 i than on r e 3 7 . As this bias 
increases, r e 2 .i deviates from the sub-pixel mean, while 
remains relatively close to the sub-pixel mean, thereby 
resulting in negative and increasing Ar ej 37 _ 2 . i . Together 
Figures 16 and 17 support the hypothesis that the plane- 
parallel r e bias is the primary reason for the systematic bias 
in the r e , 2.1 retrievals, leading to significant Ar ei 3 . 7 _ 2 .i 
values. It should be mentioned here that Marshak et al. 
[2006] found that r 6 j 2.1 retrievals tend to decrease with 
increasing horizontal scale, which is opposite to what we 
found here. The difference between these two studies might 
be attributable to the use of different LES fields, different 
levels of complexity in cloud microphysics (a constant r e = 
10 /rm for the whole cloud field was used by Marshak et al. 
[2006], while more realistic r e fields from a bin microphys- 
ics scheme are used in this study), or different configurations 
in radiative transfer simulations. 

[47] To summarize the lessons learned from the above 
case studies, we use the following equation to qualitatively 


describe the impact of various factors on MODIS r c 
retrieval: 

r e , A = r* A + A^ + A^, (2) 

where the subscript A indicates the spectral band (i.e., 2.1 /rm 
or 3.7 //m) used in the retrieval. In this equation, we use the 
r* x term to denote the “ground truth,” the A r]P x term to 
represent the impact of the 3-D radiative effects (i.e., illu- 
minating or shadowing effects) on r e retrieval and the A r p eX 
term for the plane-parallel r e bias. First, it is important to note 
that even the “ground truth” r* x is dependent on the spectral 
band owing to, for instance, the spectral difference in vertical 
weighting. Second, our case studies demonstrate that the 
A and A r pp terms and their relative impact on r eX depend 
on both spectral band and spatial resolution. At high resolu- 
tion (i.e., LES resolution), the A r^ p x term is small and Ar, ! / A 
term dominates, leading to substantial difference between 
1-D and 3-D retrievals (see Figure 13). Note that the sign of 
Ap e D x term could be either positive or negative, depending 
on the nature of the 3-D effect (i.e., illuminating or sha- 
dowing). From a spectral perspective, because the 3-D 
effects tend to impact r e21 more than r ej3 7 , the absolute 
value of ApJl. 1 tends to be greater than A 7 ^ 3. 7 (see 
Figure 14). At high resolution the sign of Ar e , 3 . 7 - 2 .i could 
be either positive or negative (see Figure 5). When cloud 
reflectance is aggregated to MODIS resolution, the cancel- 
lation of opposing 3-D effects reduces the absolute value of 
the A term. At the same time, the A r^ p x tenn from the 
plane-parallel r e bias becomes significant, especially for 
heterogeneous pixels. Note that because the sub-pixel scale 
variability of r e is significantly smaller than that of r (see 
Figure 12), the sign of the A r p x term tends to be positive 
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Figure 13. The difference between 3-D and 1-D r retrieval versus the difference between (a) 3-D r e , 2.1 
and 1-D r e2 . 1 at LES resolution, (b) 3-D r e37 and 1-D r e37 at LES resolution, (c) 3-D r e2 .\ and 1-D 
r e2 .i at MODIS resolution, and (d) 3-D r e37 and 1-D r e 3 7 at MODIS resolution. 


40 


(see Figure 1). From a spectral perspective, because the 
3.7 /jm look-up-table is more orthogonal than 2.1 /xm look- 
up-table, Ar^.i tends to be larger than A r^ 3 . 7 (see Figures 1 
and 1 1). As a result, at MODIS resolution the Ar e 3 7 _ 2 .i is 
biased toward negative values (see Figures 6 and 15) and 
generally decreases with increasing sub-pixel heterogeneity 
(see Figure 7), although for individual pixels Ar 6j3 , 7 _ 2 .i could 
be either positive or negative depending on the Ar]” term. 
Finally, it should be noted that the A and Ar(( 7 terms arise 
only when a cloud field has significant horizontal heteroge- 
neity. Therefore, these terms cannot explain why some very 
homogenous pixels in MODIS operational retrievals (i.e., 
upper-left part of Figure 2) have large negative values of 


Ar e 3 7 _ 2 j (about —5 /xm). An explanation for such pixels is 
left for future studies. 

6. Conclusions and Future Work 

[ 48 ] In this study, we develop a MODIS cloud property 
retrieval simulator based on the combination of a large-eddy 
simulation model and radiative transfer models. Using this 
simulator, MODIS retrievals of shortwave optical depth (r) 
and droplet effective radii using 2.1 or 3.7 /im radiances 
(r e . 2 .i an d r e 3 7 ) are simulated at different resolutions based 
on both 1-D and 3-D radiative transfer for several LES cases. 
The effects of drizzle and cloud horizontal inhomogeneity 
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a) SZA=50° LES resolution b) SZA=50° MODIS resolution 



Figure 14. The difference between 3-D and 1-D r retrieval versus 3-D Ar e 3 7 _ 2 .i for (a) SZA = 50° at 
LES resolution, (b) SZA = 50° at MODIS resolution, (c) SZA = 20° at LES resolution, and (d) SZA = 20° 
at MODIS resolution, based on the retrievals at nadir viewing direction from the ATEX clean case. 


on the r e retrievals in these cases are investigated. It is found 
that at high resolution (^100 m) 3-D radiative transfer 
effects, such as enhanced illumination and shadowing, affect 
r e , 2 .i stronger than r e , 3.7 probably attributable to weaker 
absorption at 2.1 /tm. As a result, the illumination effect 
tends to result in positive Ar e 3 7 _ 2 .| (= r e 3 7 — r 6 j 2 . 1 ) and the 
shadowing effect tends to result in negative Ar e 3 7 _ 2 j. 
However, because of a balance between these two opposing 
effects, r e2 .i and r c 3 7 agree reasonably well at high reso- 
lution, with no systematic bias between the two. At MODIS- 
like resolution (^800 m), r e2 . 1 is found to become 


systematically larger than r e37 and the difference is seen to 
increase with r e> 2.1 and the sub-pixel inhomogeneity index 
H a , consistent with the trends found in operational MODIS 
observations. This difference is unlikely a direct radiative 
result of drizzle here because the removal of drizzle in the 
radiative computations has little impact on the retrievals, 
including Ar e 37 _ 2 .i. It is also found that opposing 3-D 
radiative effects tend to cancel each other out at MODIS-like 
resolution, resulting in a weaker net impact of 3-D effects on 
r e retrievals. Finally, it is found within pixels of MODIS-like 
resolution that cloud r generally varies more strongly than 
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Figure 15. MODlS-resolution Ar e 3 7 _ 2 .i averaged over all 
cloudy pixels with r > 5 based on 3-D radiative transfer sim- 
ulation for the ATEX clean case at different solar and view- 
ing angles and different retrieval resolutions. 

r e . The strong sub-pixel variability in r gives rise to an 
plane-parallel r e bias that results from averaging over non- 
orthogonal r e and r retrieval curves in the 2 -channel (visible 
and SWIR) radiance phase space, and this heterogeneity bias 
is largely responsible for the systematic difference between 
r e> 3 7 and r e2 \ at MODIS-like resolution 
[49] Our results have several implications. First, it is evi- 
dent from the our analysis, as well as previous studies 
[Boeke, 2009; Hayes et al, 2010; Seethala and Hoi~vath, 
2010; Zhang and Platnick, 2011], that the cloud horizontal 
inhomogeneity effect can cause substantial errors in MODIS 
r e2 1 retrievals. The MODIS r e retrieval based on the 1.6 //in 
band (r e \ 6 ) is not considered in this study but, given that 
water absorption is even weaker in this band, we expect that 
the r e [ 6 retrieval faces the same problem. Further investi- 
gations are needed to understand the temporal and geo- 
graphical distributions of this error in the MODIS cloud 
product, and its implications for climate and aerosol indirect 
effect studies. Second, it is shown that the inhomogeneity 
index H a can be used to assess the magnitude of such errors, 
and therefore the quality of MODIS operational r c 2 . 1 
retrievals. In fact, in the upcoming MODIS Collection 6 
(C 6 ) cloud product, it has been planned for the index H a to 
be reported as part of the cloud mask product (MOD35). 
Finally, the large negative A/' e 3 7 _ 2 .i values associated with 
large H a index (i.e., points in the upper right comer of 
Figure 2) are likely the result of the cloud inhomogeneity 
effect. This finding warrants caution for the use of these 
retrievals in cloud vertical structure or drizzle detection 
algorithms based on Ar e 3 7 _ 2 1 [Chang and Li, 2002, 2003; 
Kokhanovskv and Rozanov, 2011]. 


[50] Several questions also arise from this study that are 
worthy of further investigation. For instance, is it possible to 
reduce the impact of cloud inhomogeneity on the r e2A 
retrieval? Note that the MODIS 2.1 /tm band has the native 
resolution of 500 m. For the limited cases in this study, r e 3 7 
and r e 2 ! agree better at higher resolution. Therefore, it is 
worth exploring whether the magnitude of Ar e>3 7 _ 2 [ is 
substantially reduced if retrievals are made at 500 m reso- 
lution. Such investigations would also provide useful infor- 
mation for other satellite instruments, for instance the VIIRS 
(Visible/Infrared Imager Radiometer Suite) on the Suomi 
NPP (National Polar-orbiting Partnership) mission, which 
has a resolution of about 375 m. Another intriguing question 
is whether Ar e 3 7 _ 2 3 contains any useful information for 
cloud vertical structure retrieval or drizzle detection. Note 
that a significant portion of the large negative Ar e3 7 _ 2 A 
values are found to be associated with small H a index in the 
MODIS observations (i.e., the pixels in the upper left of 
Figure 2). Such situations are not represented in our simu- 
lation cases. It seems that the cloud inhomogeneity effect is 
small for such pixels. It remains unclear and warrants further 
investigation whether these large negative Ar e 3 7 _ 2A values 
associated with small H a reflect the cloud vertical structure 
[Chang and Li, 2002, 2003; Kokhanovsky and Rozanov, 
2011] or other factors. Finally, future work is also needed 
to study the impact of cloud horizontal inhomogeneity on ice 
cloud microphysics retrievals and on other remote sensing 
techniques, for example the infrared split window method 
[Inoue, 1985; Prabhakara et al., 1988]. 




MODIS resolution x retrieval MODIS resolution r«, retrieval 




MODIS resolution x retrieval MODIS resolution r,j., retrieval 

Figure 16. Comparisons between mean retrieval values at 
LES resolution (i.e., 100 m) with MODIS-like resolution 
retrievals (i.e., 800 m) based on averaged radiance for (a) r 
and (b) r e2A using the 1-D radiative transfer simulation, 
(c, d) Same as Figure 16a and Figure 16b, respectively, but 
are based on 3-D radiative transfer simulations. 
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mean of LES resolution h n - MODIS resolution r e , u retrieval 


Figure 17. The x axis corresponds to the difference between the sub-pixel means of r e2 .i and the r e2 .\ 
retrievals at MODIS resolution. They axis corresponds to the difference Ar e 3 7 _ 2 .i at the MODIS resolu- 
tion. The figure shows that r e2 A heterogeneity bias is correlated with Ar e3 7 _ 2 A . 
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